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Using the numerical approach for a study of the thermodynamic properties of the nonuni- 
on ■ 

form one-dimensional spin- \ isotropic XY model in a transverse field we examine different 
lattice distortions to reveal which spin-Peierls phases are realized in the magnetic chain at zero 

en 

^r) ' temperature in the presence of external field. 
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An interest in the theoretical study of quantum spin chains exhibiting spin-Peierls phases has 
incredibly grown since the discovery of the first inorganic spin-Peierls compound CuGe03 (for a 



T3. 

review see 1 ). Although the quantum Heisenberg model is usually used as an appropriate model 

o 
o ■ 

to describe the spin-Peierls phase transition in the available materials some generic features can 
^ be clarified within the framework of the simpler spin-i isotropic XY chain (see 2 ' 3 and references 



therein). The latter spin model can be reformulated using the Jordan- Wigner transformation 
as a one-dimensional model of tight-binding spinless fermions. As a result, the exhaustive 
analytical and numerical analysis of different properties of the model becomes possible. 

In what follows we analyze a stability of various spin-Peierls phases at zero temperature in 
the presence of external field. For this purpose we consider a nonuniform spin- ^ isotropic XY 
chain in a transverse field defined by the Hamiltonian 

N-l N 

H=Y,Jn («+i + «+i) + n E <• (!) 

n=l n=l 

Here J n = J (1 + 5 n ) is the nonuniform exchange interaction between the sites n and n + 1 and 
the sequence of parameters 6%, . . . , <5jv-i = {$n} defines a certain lattice distortion. 0, is the 
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Figure 1: £(<5,p, ft) — £(0,p, ft) vs. 5 at different ft, a 
0.000, 0.005, ..., 0.100 (from bottom to top); (b) p = 1.9, ft 
0.100 (8); (c) p = 2.1, ft = 0.060 (1), 0.065 (2), . . . , 0.100 (9). 
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value of a uniform external field directed along z axis. Using the numerical approach described 
in detail in we calculate the ground state energy Eo({5 n } , ft) = Neo({5 n } , ft) of the spin chain 
(|l|). The ground state energy does not depend on the sign of J; we fix the units putting in what 
follows | J\ = 1. We accurately analyze the finite-size effects to be sure that our results pertain 
to the thermodynamic systems. Taking N = 1000 we no more observe the finite-size effects 
which are nicely pronounced when N = 100 or less. 

In the adiabatic treatment of the spin-Peierls instability at zero temperature one should 
examine the total energy which consists of the magnetic part Eo({5 n } ,ft) and the elastic part 
aJ2 n =i ^n f° r different lattice configurations {5 n }. Here a is the parameter which measures the 
lattice stiffness. The total energy per site will be denoted by £({8 n } ,ft). 

To begin with, we assume a = 0.5 and examine the dimerization ansatz 5 n = — 5(— l) n , 
< 5 < 1. £(5, ft) — £"(0, ft) is shown in Fig. la. Inspecting the displayed curves one concludes 
that i) if ft does not exceed ft a (~ 0.035) only the dimerized phase occurs (with the dimerization 
parameter 5* ~ 0.065); ii) if ft exceeds ft a but does not exceed ft^ (~ 0.045) both the dimerized 
and uniform phases are possible, however, the former phase is favorable; iii) if ft exceeds ft;, 



but does not exceed Q c (~ 0.065) both the dimerized and uniform phases are possible, however, 
the latter phase is favorable; iv) if £1 exceeds Q c only the uniform phase is possible. Such a 
behavior of the total energy corresponds to a scenario of the first order phase transition driven 
by the external field. Within the frames of the adopted ansatz for 8 n we are restricted to the 
dimerized and uniform phases and cannot judge about a possibility of more complicated lattice 
distortions. 

To discuss whether the total energy can be lowered by another (not uniform) lattice pattern 
as the field increases we introduce a trial distortion of the form 

S n = -5 cos I — n) , (2) 

where p is the period of modulation (e.g., p = 2 yields the dimerization ansatz). A behavior 
of £(5,p,tt) for p = 1.9 and p = 2.1 as the field increases can be seen in Figs, lb and lc, 
respectively. From the displayed plots one concludes that a long-period structure does arise, if 
£1 exceeds the value about 0.06. Thus, the dimerized phase transforms into a long-period phase 
rather than into the uniform phase while the field increases. 

Further, let us clarify whether for any small field there exists such p which yields the 
total energy lower than that for the dimerized chain. For a certain a at the fixed $7 (= 
0, 0.025, 0.05, . . .) we examine the dependence £(5,p,£l) vs. 5 seeking the minimal value 
of £ (5,p, O) for different p (in a sufficiently large region) and then compare those minimal val- 
ues. As a result, we find the value of p of the most energetically favorable lattice distortion (|2|). 
It will be denoted as p* . We repeat the search of the most favorable p at the fixed f2 for different 
a. The results obtained are collected in Fig. 2. From Fig. 2a one can learn, for example, that 
for a = 0.5 the dimerized phase is favorable at least up to f2 = 0.025 whereas at £1 = 0.05 it 
is already unstable with respect to the transformation into a long-period phase (with p about 
2.06). The important conclusion that can be drawn from Fig. 2a is that the dimerized phase for 
any a persists up to a certain characteristic field (the value of which decreases as a increases) . 
Another evident observation is that at zero field the dimerized phase brings the lowest total 
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Figure 2: p* vs. a for several values of Q. 



energy for any a (that agrees with the result proved in 5 ). 

With the increase of the field the lattice parameterized by (pf) may exhibit short-period 
phases, for example, the trimerized phase 6 for which p = 3. However, as we shall see below, a 
behavior of the trimerized phase is essentially different in comparison with that of the dimerized 
phase. Fig. 3a demonstrates that the trimerized phase may really occur at the field about 0.5. 
In the same region of fields, various long-period structures are possible (Figs. 3b, 3c). Moreover, 
for any small deviation of the field from 0.5 there exists such p for which the lattice distortion 
(|2|) gives smaller energy than for p = 3. This can be seen in Fig. 2b. Thus, contrary to 
the dimerized phase, the trimerized phase does not persist with the field varying and so it 
continuously transforms into a certain long-period phase. 

It is worth noting that using the exact analytical expression 3 for eo({5 n } ,Q) for the spin- 
2 isotropic XY chain in a transverse field of period 3 one can find that there may exist two 
lattice distortions, which preserve the chain length (Si +S2 + S3 = 0), for which £({5 n } , O) may 
have an extremum. Namely, <5i = 82 (it follows from (g) if p = 3) and 5\ = —62- We check 
numerically that in both cases £({5 n } ,0) really exhibits a minimum at 0, about 0.5 (see Fig. 
4). However, due to the instability with respect to the transformation into long-period phases 
those trimerized patterns are not so important as the dimerized configuration. 
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Figure 3: £(5,p,£l) -S(0,p,Q) vs. 5 at O about 0.5, a = 0.5. (a) p = 3, = 0.485 (1), 0.490 
(2), . . . , 0.515 (7); (b) p = 2.9, n = 0.450 (1), 0.455 (2), . . . , 0.485 (8); (c) p = 3.1, n = 0.515 
(1), 0.520 (2), ..., 0.550 (8). 



To summarize, using a numerical approach we examined (within the adiabatic approxima- 
tion) a stability of various lattice distortions of the spin- i isotropic XY chain at zero temper- 
ature in the presence of external field. We found that while the field increases the favorable at 
zero field dimerized phase persists until the field achieves a certain characteristic value, at which 
a first order transition into incommensurate phase occurs. Although in a moderate field one 
can find distortions having short periods (for example, of period 3) which bring a gain in the 
total energy, any small variation of the field leads to more energetically favorable long-period 
structures. In strong fields, the uniform lattice can be expected. Clearly, since we are restricted 
to the adopted ansatz for a lattice distortion (||) we can say for sure what lattice distortion is 
not realized rather than to point out which lattice distortion should occur. Finally, let us note 
that the elaborated numerical procedures can be easily applied to other distortion patterns, 
e.g., an array of solitons 7 . 

The authors thank J. Richter and O. Zabyrannyi for discussions. O. D. acknowledges the 
kind hospitality of the Magdeburg University in the summer of 2000 when a part of this work 
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Figure 4: £({5 n } } Q) — £ ({0}, Q) vs. 5 for two lattice distortions of period 3, i.e. 5\ = 82 = ^6, 
5 3 = -5 (a) and 8 1 = -5 2 = 8,S 3 = (b), at n about 0.5, a = 0.25. ft = 0.45 (1), 0.46 (2), 
0.55 (11). 

was carried out. 
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